Climbing the Entropy Barrier: Driving the Single- towards the Multichannel Kondo 
Effect by a Weak Coulomb Blockade of the Leads 
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We study a model proposed recently in which a small quantum dot is coupled symmetrically to 
several large quantum dots characterized by a charging energy E c . Even if E c is much smaller than 
the Kondo temperature Tk, the long-ranged interactions destabilize the single-channel Kondo effect 
and induce a flow towards a multi-channel Kondo fixed point associated with a rise of the impurity 
entropy with decreasing temperature. Such an "uphill flow" implies a negative impurity specific 
heat, in contrast to all systems with local interactions. An exact solution found for a large number 
of channels allows us to capture this physics and to predict transport properties. 
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Simple models of non- Fermi liquids, such as the multi- 
channel Kondo model, have attracted much theoretical 
interest [| , especially due to their striking properties like 
zero-bias anomalies or a finite entropy at zero tempera- 
ture, S(T = 0) = ln(g), with a non-integer g 0. So 
far, however, their experimental realization has been a 
challenging task. The advent of nanostructures that can 
be designed and tuned more easily than solid state sys- 
tems is a major step towards observing these interesting 
strong-correlation effects in the laboratory. 

The multi-channel Kondo model describes a single 
spin coupled symmetrically to several independently con- 
served conduction electron "channels". Recently Oreg 
and Goldhaber-Gordon suggested that these channels 
can be realized by attaching several large quantum dots 
serving as "leads" to a single small dot. They pointed 
out that a sufficiently large charging energy, E c > Tk, 
in those leads can supress all low-energy cotunneling 
processes between them which would otherwise mix the 
channels and destroy the multi-channel physics. The ex- 
perimental realization of such a system is a demand- 
ing task as the size of the large dots has to be cho- 
sen such that the level-splitting is sufficiently small 
and, simultaneously, the charging energy is large enough, 
A L < Tf ulti <T K <E C , where T^ ulti is the Kondo tem- 
perature of the multi-channel Kondo model and Tk refers 
to the single-channel Kondo temperature for E c = 0. 
Furthermore, considerable fine-tuning, using gate volt- 
ages, is required to guarantee that all leads couple equally 
to the spin residing on the small dot. Pustilnik et al. |j] 
have recently discussed in detail how to achieve such a 
fine-tuning by calculating the conductance for small vari- 
ations in the coupling to the various channels. The idea 
that interactions in the leads can lead to multi-channel 
physics was also put forward earlier by Coleman and 
Tsvelik 5] , and received recent additional support [(J . 

What happens if the charging energy E c of the leads is 
much smaller than the Kondo temperature, E c <C Tr-? In 
this case a single-channel Kondo resonance will develop 
upon lowering T. According to conventional wisdom, the 
Fermi-liquid fixed point of the single-channel Kondo ef- 



fect can never be destroyed by small perturbations (for 
example a weak magnetic field, B -C Tk, does not pre- 
vent the formation of the Kondo resonance). Indeed, the 
powerful "g-theorem" of Affleck and Ludwig 0, @ proves 
that no small local perturbation can destabilize such a 
zero-entropy fixed point: The impurity entropy ln(g) of 
boundary conformal field theories (to which Kondo mod- 
els belong) always decreases under renormalization group 
flow. According to this theorem, a flow from the multi- 
channel to the single channel fixed point is possible, but 
not vice versa. However, the (/-theorem does not cover 
the situation under discussion, where long-range interac- 
tions induce a Coulomb blockade in the leads. 

The purpose of this Letter is to show that in such a 
situation tiny charging energies E c <§; Tk can destroy the 
single-channel Kondo effect and stabilize a multi-channel 
fixed point. As the multi-channel system is characterized 
by a finite residual entropy [S(T = 0) = ln\/2 for two 
channels] this requires a negative impurity specific heat 
in some T range. A simple argument in favor of such a 
scenario is that the existence of the single-channel Kondo 
effect implies resonant tunneling between the leads. The 
Coulomb blockade in the leads, however, prohibits such 
a resonance. The Kondo effect cannot overcome the 
Coulomb blockade, as resonant tunneling in an energy 
window of width Tk costs a charging energy of order 
E C ((AN) 2 ) ~ E c ^f, where AN are charge fluctuations 
induced by the resonant tunneling and is the level 
spacing in the leads. These energy costs are larger than 
Tk whenever E c > Al, suggesting that for E c ^> Al a 
multi-channel Kondo effect will form. 

We consider a model where the spin on the "small" 
quantum dot (represented by Abrikosov fermions /J) 
couples to several larger dots displaying charging effects: 
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FIG. 1: Schematic plot of the multi-channel Kondo tempera- 
Note the maximum for E c ~ Tk- 



degenerate charge states in the leads. In view of greater 
generality that will be useful in the following, we sup- 
pose an arbitrary number K of interacting leads, and we 
consider an SU(N) spin in a representation which ful- 
fills the constraint £)f =1 flf a = N/2. For N = 2, the 
last term in Eq. Q describes (up to a potential scat- 
tering term) the usual exchange coupling of a spin 1/2 
to a symmetric combination of electrons from all leads, 
ct _ = -= V5 i ct„„. Therefore the conventional one- 
channel Kondo effect develops for E c — 0. Note the sym- 
metric coupling of the spin to all leads; as mentioned 
above, this requires fine-tuning in experiments. 

In order to make contact with previous results 0-13 ■ we 
first use the perturbative renormalization group (RG) to 
discuss the case of large charging energy Tk <C E c <C D, 
where D is an ultra-violet cutoff set by the bandwidth or 
the charging energies in the small dot. For cutoffs large 
compared to E c , the Coulomb blockade does not modify 
the RG flow for the dimensionless coupling j = JNf 
(Nf is the density of states) . For a running cutoff A S> 
E c , one therefore finds within one-loop RG dj/dlnA = 
-f 



2 and at the scale A = E c the running j(A) takes the 



value j(E c ) = l/\n[E c /T K ] with T K = D, 
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j(E c ) < 1 as £ c > Tk- For A <C E c the charge on each 
lead is conserved separately and all processes which mix 
channels arc frozen out. Therefore one expects a flow 
to the /^-channel Kondo fixed point with dj/dln A = 
—j 2 /K with the flow starting at j(E c ). From this we 
can read off the multi-channel Kondo temperature 9] for 

t k <s:e c <s:d 
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For E c > _D, the multichannel equation determines the 
RG flow alone as channel mixing is suppressed, which 
leads to the usual result T£? ulti = De~ K / J . Fig. [T] shows 
the enhancement of Tg? ultl with decreasing E c . 

To analyze the physics at small E r , we introduce a 
phase (or slave rotor) representation |1CX fll| of the charg- 
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where we have set c krTa =a\ rTa exp(i9 a ) and L a ——id/d8 a . 
This effective matrix of Kondo couplings J aa > {0) = 
J exp(i9 a — i6 a i) allows to reinterpret the original idea 
of Oreg and Goldhaber-Gordon. At large E c the phases 
6 a fluctuate wildly in J aa >(0) for a ^ a' ', which leads to 
a flow towards a diagonal coupling at low energy. In this 
case, channel number is conserved and a multi-channel 
Kondo effect will develop. However, for E c equal to zero 
strictly, the phases are locked and drop from the Kondo 
coupling in Eq. J3J) , leaving a model displaying Fermi liq- 
uid properties and complete screening of the spin. 

To study the effects of a small E c -C Tk within per- 
turbation theory, we expand Eq. @ around constant 8 a . 
This expansion is valid as long as {{0 a (r) — # q (t')) 2 ) <C 1, 
which allows us to reliably calculate the correction AF to 
the free energy for T 3> E c . After some straightforward 
but tedious manipulations we obtained 
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where T(t) is the T-matrix of the symmetric channel 



and G (iu! n ) = £\ l/(i 



e k ) the bare local Green's 



function of the conduction electrons. For the last equality 
in |0J we used that wNflmT(u) = 1 for to, T < T K which 
follows from Friedel's sum rule. It is therefore valid for 
E c T Tk ■ Up to (known) prefactors of order 1 , the 
impurity entropy in this regime is thus given by 



S(T)/N « T/T K + E c /T, 



(5) 



where the first term arises from the usual Kondo effect 
and the second from the charge fluctuations described by 
Eq. Q . This proves that the impurity entropy will show 
a minimum at the scale Tg = \JE c Tk and implies a nega- 
tive impurity specific heat for E c <C T <C T$ (the specific 
heat of the total system ~ T/Al remains positive as we 
assumed a neglegible level spacing -c E c ). Moreover, 
the entropy per spin reaches values of order 1 for T ~ E c 
where the expansion breaks down, opening the possibil- 
ity of a flow towards a multi-channel fixed point. When 
we take into account that the crossover at E c takes place 
deep in the strong-coupling regime, this suggests that 
the corresponding multi-channel Kondo temperature is 
directly given by E c , as shown schematically in Fig.Q] 

To be able to describe the crossover at E c and the 
physics for T < E c <C Tk we need a non-perturbative 
method. We have found an exactly solvable limit of the 
Hamiltonian Q that confirms the previous calculations, 
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and also offers a direct computational tool to describe the 
crossover from single-channel to multi-channel physics as 
a function of T and charging energy. The idea is to 
solve the problem by taking both N and K to be large 
(we recall that we have considered a generalized model 
with an SU (N) quantum spin in the dot and K = 7 AT 
interacting leads coupled to it). The technical step is 
to notice that the Kondo interaction in © can be de- 
coupled using a single bosonic field B(t) conjugate to 
E/cctq fa a kaa exp(-i0 a ). Integrating out the leads, we 
obtain the following action in imaginary time: 



f KB^B 



£ 



(d T 9 a ) 2 
4E C 



£4(8r-M)/* 



N 



dr/dr' G ° {T ^ YiVlBe-v-W.BW-U (6) 
o Jo n 



where fx is a complex Lagrange multiplier used to en- 
force the constraint on the spin size. We can introduce 
[T^| two fields Q(t,t') and Q(t,t') to decouple fermions 
from bosons in the last term of ©■ Then the /J and 9 a 
variables are integrated out to obtain an effective action 
S[B,Q,Q, n] which is proportional to N and therefore 
solved by a saddle-point when N — > oo. Using time- 
translational invariance and particle-hole symmetry (so 
that (i = 0), we obtain the self-consistent equations 
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where the condensate B is determined from Eq. Ijlll) and 
the correlator Gx{t) is computed from the action 
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In principle, the model i|12fl can be solved very efficiently 
by Monte Carlo , but we have chosen to simplify the 
numerics by using the spherical limit as an approximation 
to this rotor model. This is done by introducing a further 
large M expansion, where we generalize e* e = X with the 
constraint IX | 2 = 1 to M fields X t with Y^Li \M 2 = M 
(see Ref. |Tj| for details). We then obtain 

G x \iv n ) = vl/{AE c ) + \ + Q{iv n ) (13) 
Gx(t = 0) = 1. (14) 

The parameter A is determined from Eq. H14JI. reflecting 

the constraint \e l6 \ =1 in average. 

We first investigate the general properties of the sys- 
tem of equations l|7lll|) and I|13I14|I . For E c — we re- 
cover nicely the usual large- N limit of the single-channel 
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FIG. 2: Entropy S(T) for T K /E C = 400, 80, 40, 20, 8, 4 
(bottom to top) with JNf = 1 /4, 7 = 1 and a flat density of 
states of half- width D — 2 (giving Tk ~ 0.04). The minimum 
of S is located at T s « Q.T^/eITk (crosses). For T < T s 
and _B C Tk, the entropy is a function of T/E c only (see 
inset), which shows that E c can be identified with the multi- 
channel Kondo temperature T^ ultl in this limit. The dotted 
line in the inset denotes the exact value for S(T = 0) taken 
from Ref. Note that in the regime 0AE c < T < Ts, the 
impurity specific heat C = TdS/dT is negative. 



Kondo model [T^ , since Gx (t) = 1 follows from Eqs. I|13t - 
I14|l in this case. For E c ^ 0, this behavior also holds at 
E c <C T, reflecting the fact that the single-channel Kondo 
fixed point controls the regime E c -C T -C Trt- However, 
at T = a low-frequency analysis of our integral equa- 
tions |lll Il2| shows the appearance of universal power 
laws characteristic of the multi-channel fixed point |12| . 
G f {iu) ~ ^Iuj] 1 /^ and G x - l/|w|V(i-Hy). 

Our equations allow also to calculate physical quan- 
tities at intermediate coupling. Fig. [21 shows the im- 
purity entropy, i.e. the difference of the total entropy 
and the entropy in the absence of the spin. As pre- 
dicted by Eq. JSJ, there is a minimum in S(T) at the 
scale Ts — \/E c Tk and the specific heat is negative for 
E c < T < Ts- Close to the low-T fixed point the specific 
heat is positive. Accordingly, the entropy drops below E c 
and for T — > reaches the finite value characteristic for 
the if-channel fixed point The scaling plot in the 

inset of Fig. [21 shows that E c can be identified with the 
multi-channel Kondo temperature Tj^ ultI if E c <C Tk, as 
discussed above. For E c > Tk, we find only the free-spin 
solution with B = and no multi-channel physics. This 
is consistent with Eq. JSJ , which shows that T|? ultI drops 
rapidly to for K — > 00 if E c > Tk- For finite K and 
low T , multi-channel behavior is of course maintained for 
arbitrarily large E c . 

Although this thermodynamic analysis provides in- 
teresting insights into the model considered here and 
the general question of how long-range interactions can 
destabilize zero-entropy fixed points, it is almost impos- 
sible to measure the entropy of quantum dots. Therefore 
we have also calculated the conductance through the dot 
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FIG. 3: Conductance G(T) for E c = 0.0001, 0.0005, 0.001, 
0.002, 0.005, 0.01 (top to bottom) as in Fig. H For T -> 
and 7 = 1, G approaches the universal value G = nGo/4. In 
contrast to the entropy (inset of Fig. |5J, large non-universal 
corrections spoil the one-parameter scaling with T/E c . Inset: 
Experimental setup suggested in Ref. Q- Two leads (1l and 
1r) are connected to a voltage source, the other "leads" 2, 3 
are large quantum dots with a charging energy E c and small 
level spacing Al <§C E c , prohibiting charge transport between 
1 and 2, 3. Therefore, a mapping onto our Hamiltonian Q 
with K = 3 is possible, where channel 1 arises from the even 
combination of electrons 3, 4] in 1l and 1r. However, fine- 
tuning is required Q to obtain symmetric coupling to all leads 
at lowest energies. 

using a setup suggested by Oreg and Goldhaber-Gordon 
which is sketched in Fig. EJ In a situation where the 
odd combination of electrons in leads 1 l and 1 decou- 
ples from the dot, the linear conductance can be calcu- 
lated from the imaginary part of the T-matrix 0, : 

° = 2^ / d-^ImGoMlmTH (15) 
T(r) = B 2 G x {T)G f {r). (16) 

For E c <C T <C Tk, the single-channel Kondo effect 
results in resonant scattering among the K equivalent 



leads and therefore in a conductance G « (N/K)Gq = 
(1/7)G with G = e 2 /{2nh). For T < E c « T™ ulti , the 
conductance is governed by the multi-channel T-matrix 
which has been calculated from conformal field theory 
by Parcollet et al. ^3 for arbitrary K and N. For 
K = N = 2 one obtains G — Go, while in the large 
N limit one gets G = G tt/(2 + 27) tan[7r/(2 + 27)]. 

We believe that all the results obtained in the previ- 
ous large N and K limit are qualitatively valid for the 
experimentally relevant case N = 2 and K — 2 or 3. 
First, this is indicated by our perturbative expansion Q 
for T 3> E c and -E c <C Tk which proves the existence of 
the minimum in S(T) for arbitrary values of N and K . 
Second, the Coulomb blockade of the leads makes the 
multi-channel fixed point obviously stable against inter- 
lead tunneling. Finally, we have verified our scenario for 
K = 2 in a strong-coupling expansion of the Hamilto- 
nian m, taking first the limit J — > 00 (and therefore 
Tk — > 00) and analyzing the resulting model for large 
E c -c J- In this limit we recover a 2-channel Ander- 
son model which can be mapped via a Schrieffer- Wolff 
transformation to the two-channel Kondo model. 

In conclusion, we have shown that tiny charging ener- 
gies in the leads, E c <C Tk, can destabilize the single- 
channel Kondo effect and induce a flow towards the 
multi-channel Kondo fixed point. While in systems with 
local interactions the g-theorem 0, guarantees that 
the impurity entropy always decreases with decreasing 
T, in our case it will rise for E C <T< \/E c Tk- Our ob- 
servation that the multi-channel Kondo fixed point can 
be stabilized even by small charging energies should help 
to realize the experimental setup proposed by Oreg and 
Goldhaber-Gordon Q. To obtain a high multi-channel 
Kondo temperature TjJ? ultI , parameters with E c sa Tk 
seem to be the most promising candidate (see Fig.0. 
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